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Abstract 

In an effort to validate computational fluid dynamics procedures for calculating 
flows around circulation control airfoils, the commercial flow solver FLUENT was uti- 
lized to study the flow around a general aviation circulation control airfoil. The results 
were compared to experimental and computational fluid dynamics results conducted at 
the NASA Langley Research Center. The current effort was conducted in three stages: 
(i) a comparison of the results for free-air conditions to those from experiments, (ii) a 
study of wind-tunnel wall effects, and (iii) a study of the stagnation-point behavior. 
In general the trends in the results from the current work agree well with those from 
experiment, some differences in magnitude are present between computations and ex- 
periment. For the cases examined, FLUENT computations showed no noticeable effect 
on the results due to the presence of wind-tunnel walls. The study also shows that the 
leading-edge stagnation point moves in a systematic manner with changes to the jet 
blowing coefficient and angle of attack, indicating that this location can be sensed for 
use in closed-loop control of such airfoil flows. 


Nomenclature 

A area 

b wing span 

c chord 

Cd drag coefficient 

Ci lift coefficient 

C m pitching moment coefficient about quarter chord 
Cfx momentum coefficient 
h slot height 

m mass flow rate 

M Mach number 

P pressure 

q dynamic pressure 
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R gas constant for air 

r radius of coanda surface 

Re Reynolds number 

T temperature 

U velocity magnitude 

w slot width 

a angle of attack 

7 ratio of specific heats 

p viscosity 

p density 

Subscripts 

duct stagnation conditions inside plenum 
fc conditions at flow-control boundary 

00 freestream conditions 

J slot-exit conditions 

1 Introduction 

Recent research in the Applied Aerodynamics Group at NCSU has led to the devel- 
opment of an automated cruise-flap system. 1,2 The cruise flap, introduced by Pfen- 
ninger, 3,4 is a small trailing-edge flap which can be used to increase the size of the 
low-drag range of natural-laminar-flow (NLF) airfoils. The automation is achieved 
by indirectly sensing the leading-edge stagnation-point location using surface pressure 
measurements and deflecting the flap so that the stagnation-point location is main- 
tained at the optimum location near the leading-edge of the airfoil. Maintaining the 
stagnation point at the optimum location results in favorable pressure gradients on 
both the upper and lower surfaces of the airfoil. With such a cruise-flap system, the 
airfoil is automatically adapted for a wide speed range. This automated cruise-flap 
system was successfully tested in the subsonic wind tunnel at NCSU. 2 

While the use of a cruise flap on an NLF airfoil results in low drag over a large range 
of speeds, there is a need for a revolutionary approach that integrates the achievement 
of significantly lower drag over a large range of operating speeds with the capability 
for generating very high lift at takeoff and landing conditions. Toward this objective, 
it is of interest to study an approach that integrates aerodynamic adaptation with 
the well-established high-lift capability of circulation control (CC) aerodynamics. This 
aerodynamic adaptation carries with it the possibility for significant skin-friction drag 
reductions through extensive laminar flow in addition to the high-lift benefits of CC 
aerodynamics. Figure 1 illustrates the overall concept. In a manner similar to a cruise 
flap, it is believed that by utilizing this stagnation-point sensing scheme, an adaptive 
CC airfoil can achieve extensive laminar flow over a large lift-coefficient range. 

As a first step toward the long-term goal of studying an adaptive CC airfoil, the 
current effort was undertaken for establishing and validating computational fluid dy- 
namics (CFD) analysis procedures for blown-trailing-edge airfoils. The CFD package 
used for this work was the FLUENT flow solver. The results are compared to CFD 
and experimental data obtained from a recent study by Jones et al. 5 of a General 
Aviation CC (GACC) airfoil conducted at the NASA Langley Research Center. Since 
previous CFD studies on this airfoil did not include tunnel walls, the current CFD 
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study also includes an investigation of the effect of tunnel walls on the solution. In 
order to provide a foundation for the adaptive CC airfoil concept, the effects of CC on 
the leading-edge stagnation-point location was also examined in the current work. 

2 Approach 

2.1 CFD study 

The commercial flow-solver code FLUENT version 6.1 was used in the current research. 
Grid generation was performed using GAMBIT, which is the preprocessor packaged 
with the FLUENT code. These codes were used to study two separate cases. The 
first case involves the examination of the GACC airfoil in free air with the objective of 
comparing the FLUENT results to CFD and wind-tunnel results presented in Ref. 5. 
The second case involves simulations of the GACC airfoil in the Basic Aerodynamic 
Research Tunnel (BART) to examine the influence of tunnel walls on this particular 
airfoil. Results from FLUENT were obtained for a matrix of 15 data points for both 
the cases. 

2.1.1 Geometry and grid details 

The geometry chosen for the current research was the General Aviation Circulation 
Control (GACC) airfoil, designed by Jones. 5 The GACC airfoil was derived from a 
17% GAW(l) airfoil by modifying the trailing edge to incorporate a 2% r/c coanda 
surface and is shown in Fig. 2. 

For the first study, a circular computational domain (Fig. 3) was generated that 
extends to approximately 20 chord lengths in all directions and is comprised of 132,762 
cells. For the study of wall effects, a second grid was generated to include the wind- 
tunnel geometry and is shown in Fig. 4. The experiments were conducted by Jones et 
al. 5 in the BART wind-tunnel which is located at the NASA Langley Research Center 
in Hampton, Va. The BART tunnel has a physical test-section size of 28” x 40” x 120”. 
The GACC model chord length was 9.4” with angle of attack changes made about the 
half-chord location. The details of the experimental setup are given in Ref. 6. For the 
computation with walls, a separate grid was generated for each angle of attack, each 
of which is comprised of 123,602 cells and extends to 20 chord lengths upstream and 
downstream of the airfoil. 

The grids for all of the analyses are hybrid unstructured grids. The domains consist 
of an unstructured grid far from the airfoil in order to reduce the number of cells 
and structured grid near the airfoil in order to maintain good resolution through the 
boundary and shear layers. 

2.1.2 Solver settings 

For the current study, the steady, coupled, and implicit solver settings with node-based 
discretization scheme were selected. The coupled solver was chosen for two reasons. 
First compressibility effects need to be modeled, as the Mach number at the slot exit 
can often approach the sonic condition as the blowing rate is increased. Secondly, 
the FUN2D code has a compressible solver and because the results from the current 
study were compared with FUN2D results, a compressible solver was also used for the 
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FLUENT analysis. There was an attempt to run these problems with the segregated 
(decoupled) solver using very low relaxation factors, however it was found that for 
the cases with larger blowing rates the solution began to exhibit an unsteady effect 
after a few thousand iterations. In order to compare with the FUN2D 5 results, the 
one-equation Spalart-Allmaras turbulence model was chosen. 

2.1.3 Boundary conditions 

FLUENT does not allow the user to input freestream Mach number and Reynolds 
number directly. Instead, the freestream velocity and operating pressure were calcu- 
lated using Eqs. 1-3 and provided as inputs for the analyses. The Mach and Reynolds 
numbers were set to 0.1 and 533,000, respectively, to match those used in Ref. 5 


r oo = M 00 ^ 1 RT 00 

(1) 

Refioc 

l> °°~ U x c 

(2) 

Poo = PooRToo 

(3) 


An approximate method was developed to estimate the required velocity at the 
flow control boundary (C// c ) to achieve a desired C^, C \i desired - This method assumes 
incompressible flow throughout the duct, and was derived by solving the continuity 
equation. The equation for U f c from this approximate method is given in Eq. 4. 

Ufc = Uoo \l (4) 

Once FLUENT converged, an integration was performed across the slot exit as 
shown in Eq. 5 to obtain the actual of the jet at the slot. This T' /( , however, 
is different from Cfj, deeired because the U f c for the latter is set using an approximate 
method. 


I slot pV dy 

^integrated ~ 1 ^ > 

Furthermore, in order to be consistent with the methods used for calculating 
in Ref. 5, all of the values presented in this paper were calculated using isentropic 
flow relations. 5 The equations for this procedure are given in Eqs. 6-8. In order 
to determine how close the isentropic C\ L is to the integrated C^, the two values are 
compared in Fig. 5 for several cases. The values indicated along the x-axis are values 
calculated using the isentropic relations. Values for C M on the y-axis were computed by 
integrating the flow across the slot exit. The solid line in Fig. 5 indicates where the data 
points would lie if the two methods generated the same values for . The symbols are 
representative of the actual values calculated using FLUENT and isentropic relations. 
Although the differences are very small, care must be taken to ensure consistency in 
the CFD solutions and experiments. 


m = pjUjAj 
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3 Results 


The results from FLUENT predictions for the GACC airfoil are presented in three 
parts. In the first part, the prediction for the GACC airfoil in free-air conditions are 
compared with the results presented in Ref. 5. In the second part, the predicted results 
for the GACC airfoil with tunnel walls are presented and compared with the free-air 
results. In the third part, the effect of a and on the leading-edge stagnation- point 
location are presented and discussed. 

3.1 Results for free- air conditions 

In this part of the study, FLUENT results for free-air conditions are compared with 
CFD and experimental results from Ref. 5. The comparison is illustrated using Ci-a 
curves in Fig. 6. The results from FLUENT analyses consist of a matrix of 15 data 
points for a =— 5, 0, and 5 deg and = 0, 0.008, 0.024, 0.047, and 0.078 and are 
presented in Fig. 6 using red dashed lines and square markers. The wind-tunnel results 
from Ref. 5 are presented as blue markers with best-fit lines in Fig. 6 for several angles 
of attack and for = 0, 0.007, 0.015, 0.025, 0.041, and 0.060. The values of C \i for 
the FLUENT results differ from those for the results of Ref. 5 because of the difference 
between the actual and the desired Cp when using the approximate method in Eq. 4 
for estimating the Uf c using incompressible-flow equations. 

Although the values of for the FLUENT results do not match those for the 
results of Ref. 5, it is clear that the trends and most of the predictions for the Q are 
close to those from Ref. 5. In particular, the FLUENT predictions for = 0, 0.008, 
and 0.047 agree quite well with the results for similar values of from Ref. 5. Two 
discrepancies between the FLUENT predictions and those from Ref. 5 are apparent: 
(i) for the C ^ = 0.024 and (ii) for C = 0.078. The reason for the first discrepancy 
in the results is attributed to the incorrect prediction of the jet-separation location on 
the Coanda surface for = 0.024. The apparent discrepancy in the results for = 
0.078 is attributed to nonlinear effects at the high blowing rates and the fact that the 
highest blowing rate in the results of Ref. 5 is for C M = 0.060. 

The flow-field data for the FLUENT results are presented in two separate parts. In 
the first part, the effects of increasing C /( for a constant angle of attack is presented. 
The second part examines the effect of angle-of-attack changes and their influence on 
the CC airfoil for a constant C M . The flow-field data is presented as pressure contours 
and streamline plots; these aid in the understanding of the effects of CC on the flow 
over the airfoil. 

The first part of the flow-field data is shown in Figs. 7(a)-(c). It can be seen that 
as the blowing rate is increased the streamlines become more curved — an indication 
of increased circulation. The second part of the flow-field data is shown in Figs. 8(a)- 
(c) and Figs. 9(a)-(c) to illustrate the effects of changing the angle of attack while 


817 


holding blowing rates constant. The results are presented for two blowing rates: the 
mild blowing case C M = 0.047 and the highest blowing rate C ^ = 0.078. The results 
show that changes to have a significant effect on the jet-separation location and 
the resulting C/. In comparison, changes to a have a much smaller effect on the jet- 
separation location. 

3.2 Wind-tunnel wall effects 

In this sub-section, the FLUENT results for the GACC airfoil with the effect of wind- 
tunnel walls are presented. Figures 10-12 show the influence of the wall on the CFD 
solution. These figures present the predicted Ci as a function of for a = 0, 5, and 
—5 deg respectively. Figure 10 also includes a comparison to results for the FUN2D 
study 5 for a = 0 deg, the only angle of attack for which the FUN2D results were 
presented in Ref. 5. Figures 10-12 indicate that the presence of walls has very little 
influence on the CFD solution. The solution, including tunnel walls, consistently show 
that for low blowing coefficients, the Ci values are predicted to be lower than those 
without walls. However, at the largest blowing coefficients, the trend reverses and Ci 
values with walls are predicted to be higher than those without walls. 

3.3 Stagnation-point location 

The motivation for examining the stagnation-point behavior is that the stagnation- 
point location was used successfully in earlier research 1,2 for closed-loop control of a 
trailing-edge flap. It was, therefore, desirable to examine the CFD solutions for the 
CC airfoils to see if there was any evidence that would suggest that a similar approach 
could be extended for use with CC airfoils. 

Stagnation-point location, measured as an arc length from the jet exit around the 
upper surface of the airfoil, as a function of Ci is presented in Fig. 13. Each line 
in Fig. 13 represents a different blowing rate and for each blowing coefficient there 
are three points that correspond to three different angles of attack(— 5,0, and 5 deg). 
From Fig. 13 it can be seen that the stagnation point moves in a predictable manner, 
both with angle of attack and with changing blowing rate. This behavior provides 
an indication that the stagnation-point location can be used as a means to develop 
closed-loop control of the jet on CC airfoils. 


4 Conclusions 

The results from a two-part CFD study using the FLUENT flow solver have been 
presented. Results of the first study show that while the FLUENT predictions do 
not match the CFD and experimental results of Ref. 5 exactly, the overall trends 
are followed very closely. Throughout the range of blowing coefficients, FLUENT 
consistently predicted a slightly lower overall lift coefficient. 

In addition, a study was performed on the influence of wind tunnel walls on the 
CFD solution. For low blowing coefficients, it was found that the lift is predicted 
to be lower for the cases with walls. The trends are reversed for the higher blowing 
coefficients, for which the cases with walls yield a higher predicted lift. Although the 
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solutions are different, the differences are small, and could as well be attributed to 
differences in the grids rather than the actual presence of walls. 

The influence of circulation control on the leading-edge stagnation point location 
was examined. It was shown that changes in blowing rate and angle of attack result in 
systematic changes to the stagnation-point location. This observation indicates that it 
is possible to use a closed-loop control system by sensing the stagnation-point location. 
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Figure 1: Illustration of the NCSU concept of an adaptive circulation control airfoil. 



Figure 2: General Aviation Circulation Control (GACC) airfoil geometry used in the current 
research. 



Figure 3: Grid generated for FLUENT comparison to FUN2D. 
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Figure 4: Grid generated for FLUENT study of wall effects. 



Figure 5: Comparison of C Minte9rote<i with C\ lisentropic for a - 0; the straight line is included to 
indicate deviation from a perfect correlation. 
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C p : -0.77 -0.53 -0.30 -0.06 0.17 0.41 0.64 0.88 




(a) - 0.000 


>• 


>- 


(b) - 0.047 



(c) Cf, = 0.078 

Figure 7: CC effects on the flow field at a = 0 for various values of C \ 1 . 
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- 0.77 - 0.53 - 0.30 - 0.06 0.17 0.41 0.64 0.88 



(a) a = — 5 deg 



(b) a = 0 deg 



0 v 1 0.7 0.8 0.9 1 1.1 1.2 

x X 

(c) a — 5 deg 

Figure 8: CC effects on flow field at C\ L = 0.047 for various values of a. 
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C p : -0.77 -0.53 -0.30 -0.06 0.17 0.41 0.64 0.88 



(a) a — —5 deg 



(b) a — 0 deg 



(c) a — 5 deg 


Figure 9: CC effects on flow field at — 0.078 for various values of a. 
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Figure 10: FLUENT prediction of wind-tunnel wall effects for varying values of at a 
0 deg. 



Figure 11: FLUENT prediction of wind-tunnel wall effects for varying values of C\ t at a 
5 deg. 
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Figure 12: FLUENT 
—5 deg. 


prediction of wind-tunnel wall effects for varying values of at a 
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Figure 13: CC effects on stagnation-point location. 
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Background & Motivation 

• A recent NCSU Applied Aero Group (APA) 
accomplishment in adaptive aerodynamics 

• NCSU APA long-term goal in CC 

• Current CC project in collaboration with Dr. 
Greg Jones, NASA LaRC 




Recent NCSU APA Accomplishment 


• Aerodynamics for adaptive aircraft 


• Recent successful work on auto-adaptive 
airfoil (TE flap) 


Example: 

NASA NLF(1)-0215F at Re 
= 6 million 



Flap = -10 deg 
Flap = 0 deg 
Flap = 10 deg 
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Long-term NCSU APA Goal in CC 


• Adaptive jet-flap/CC concept 


high-speeo 





Current CC Research Project 


• In collaboration with Dr. Greg Jones (LaRC) 

• Use Fluent code to analyze the GACC airfoil 

• Compare results with LaRC experiments and 
FUN2D computations 

• Build CFD expertise in CC for follow-on 
adaptive jet flap research 
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Geometry 


• General Aviation Circulation 
Control (GACC) Airfoil 
(Supplied by Dr. Greg Jones) 

• Designed for upper and 
lower blowing 

• Conditions for landing/take- 
off conditions 

- Re M = 533,000 

- Ho = o.i 




Machanical & 


Grid Details 


• Generated using Gambit 

• Hybrid unstructured grid 

• Structured near airfoil 

• Unstructured in far-field 

• ~ 20 chord lengths 

• 132762 elements 
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Solver Details 

• Fluent version 6.1 

• Steady 

• Coupled-Implicit 

• Node-based discretization 

• Spalart-Allmaras turbulence model 



Results 

• Current Fluent work compared to Jones et 
al. (2002) 

- Experiment 

- FUN2D CFD analysis 

• Parametric study of a and Cp variations 
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Comparison of Lift Characteristics 




Parametric Studies 

• Increasing C at a = 0° 

H' 

• Increasing a at C = 0.047 

M' 

• Increasing a at C =0.078 
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Trends in Lift and Pressure Coefficient 
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Airfoil Surface Pressure Coefficient 
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Airfoil Surface Pressure Coefficient 


C = 0.078 





CC Effect on LE Stagnation Point 
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Conclusions 


• CFD computations performed for GACC with 
Fluent 


• Trends and much of data compare well with 
LaRC experiments and FUN2D 

• Typical run time of 1 .5-2.5 days per case 
(Pentium Xeon 3.0 GHz) 


• Systematic movement of LE stagnation point 
can be used for sensing and closed-loop control 




Provides foundation for adaptive jet flap/CC 
research 
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Mach contours for increasing a 
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Mach contours for increasing a 
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Typical Residual Plot 
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Cruise Flaps for NLF Airfoils 

• Trailing-edge flaps or “cruise flaps” can be used to 
shift the drag bucket of an airfoil. 

• Originally conceived by Pfenninger (1947) 


Example: 

NASA NLF(1)-0215F at Re 
= 6 million 
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Stagnation-point Sensing 

Cruise flaps work by maintaining the stagnation 
point within a small desirable region near the 
leading edge (well known) 



Example: 

NASA NLF(1)-0215F 



Sensing can be done using surface hot-film array 
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